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We present the first analytical superposition of a charged black hole with an annular disk of 
extremal dust. In order to obtain the solutions, we first solve the Einstein-Maxwell field equations 
for sources that represent disk-like configurations of matter in confomastatic spacetimes by assuming 
a functional dependence among the metric function, the electric potential and an auxiliary function, 
which is taken as a solution of the Laplace equation. We then employ the Lord Kelvin Inversion 
Method applied to models of finite extension in order to obtain annular disks. The structures 
obtained extend to infinity, but their total masses are finite and all the energy conditions are 
satisfied. Finally, we observe that the extremal Reissner-Nordstrom black hole can be embedded 
into the center of the disks by adding a boundary term in the inversion. 

PACS numbers: 04.20.Jb 04.40.Nr 98.62.Js 



I. INTRODUCTION 



Exact solutions of Einstein equations with axial symmetry play an important role in astrophysical applications of 
general relativity. In particular, disk-like configurations of matter are of great interest, since these structures are 
present in several systems, such as galaxies, nebulae and accretion disks around compact objects. In fact, there are 
still important unanswered questions in this field, namely the formation of Active Galactic Nuclei (AGN), X-ray 
transients and gamma-ray bursts, which are commonly associated with the accretion of matter onto black holes. 

It is commonly thought that these sources are able to form collimated, ultrarelativistic flows (relativistic jets). 
The exact mechanisms to explain the production of jets are still uncertain, but they probably involve the interaction 
between a rotating black hole, the accretion disk and electromagnetic fields in strong gravitational fields (see, for 
example, [H-HI and references therein). Thus, an accurate general relativistic model of an AGN would require an 
exact solution of Einstein-Maxwell field equations that describes a nonlinear superposition of a Kerr black hole with 
a stationary disk coupled to electromagnetic fields. Not even an exact solution of a stationary black hole-disk system 
has been found yet. 

Hence, the study of systems composed by a thin disk surrounding a central black hole is of great relevance in 
astrophysics and it is considered an open problem in general relativity. A lot of work has been done in the last years 
in order to obtain exact solutions of Einstein's equations corresponding to thin disk-like sources with a central black 
hole as well as to obtain a better understanding of the different aspects involved in the dynamics of these systems 
(see d, [|[ for thoroughly reviews on the subject). 
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Solutions for static thin disks without radial pressure were first studied by Bonnor and Sackfiel Q and Morgan 
and Morg an Q and with radial pressure by Morgan and Morgan Q . Other classes of static thin disks were obtained 
by P4l4j| and stationary thin disks were studied in fl5l - [T8l | . Models of thin disks with halos made of perfect fluid 
were considered in 1191 . The generalization of the displace, cut and reflect method for constructing static thick disks 
was considered by |20| . Thin disks have been discussed as sources for Kerr-Newman field [U [22j , magnetostatic 
axisymmetric fields |23l | and conformastatic and conformastationary metrics [2"H [25| , while models of electrovacuum 



static countcrrotating dust disks were presented in [26|. Charged perfect fluid disks were also studied in [27j, and 
charged perfect fluid disks as sources of static and Taub-NUT-type spacetimes in [28|, [29[ . 

However, in the context of a superposition with a black hole, if we consider disks that extend up to the event 
horizon, the matter located near the black hole will have superluminal velocities, as was shown by Lemos and Letelier 
[30M32I ] . So, in order to prevent the appearance of tachyonic matter, the thin disks must have a central hole with a 
radius larger than the photonic radius of the black hole. Then, the boundary value problem is mathematically more 
complicated and thus only very few exact solutions have been obtained. 

This kind of solutions were first studied by Lemos and Letelier [3l| for static axisymmetric spacetimes to obtain 
exact disk solutions around a Schwarzschild-type black hole. Since the static axisymmetric Einstein equations are 
equivalent to the Laplace equation and a nonlinear quadrature, this supcrsposition of solutions is possible there. The 
main properties of this kind of superposition were extensively analyzed in a series of papers by Semerak, Zacek and 
Zellerin [3314391 ] . by using numerical computation when needed. Besides these models, only two other solutions for 
static thin disks with an inner edge have been obtained, a first one with inverted isochrone disks (4fj| and a second 
one for disks with a power-law density [4l|. Also, Zellerin and Semerak (42j, by using the Belinskii-Zakahrov inverse- 
scattering method, found a stationary metric that reduces to the superposition of a disk and a black hole in the static 
limit and thus may represent a stationary disk-black hole system, but the analysis of their solution is complicated by 
the fact that the metric functions cannot be analytically computed. Furthermore, this solution involves an unphysical 
supporting surface between the black- hole horizon and the disk (43[. Finally, general class of stationary solutions 
was presented by Klein [3,|4^], using the Riemann-surface techniques, in which physically acceptable black hole disk 
systems can be found. At this point it is worth to mention that a common feature of all the above mentioned solutions 
is that their metric functions can not be fully analytically computed and most of them present singularities at the 
inner edge of the disk. Thus, the analysis of their physical and mathematical properties is very complicated. Recently, 
the first fully integrated exact solution of the Einstein eq uations corresponding to the superposition of an annular 
dust disk with a neutral black hole was founded in [461 147| . overcoming most of the problems mentioned above. 

In this paper we present the first superposition of a charged black hole with an annular disk made of extremal 
dust. Although one may intuitively expect that astrophysical objects are neutral, the possibility that compact objects 
could actually contain a non- vanishing net charge was first pointed out in [48l |43 | , where the authors modeled a star 
as a gas of positive ions and electrons and concluded that, due to their greater kinetic energy, the electrons tend 
to escape from the star more often than the ions. The star will then acquire a net positive charge. The process 
will be carried on until the electric field induced in the star stops more electrons from escaping. Recently other 
mechanisms to induce electric charge into compact objects, in particular into black holes, have been proposed (see 
for example (Hoj). In order to obtain the superposition, we first construct two infinite families of solutions of the 
Einstein-Maxwell equations in conformastatic spacetimes that represent disk-like configurations of matter with a 
central hole and charge density equal to the mass density, in such a way that the electric repulsions and gravitational 
attractions are in exact balance. This kind of equilibrium configuration of matter, has been called by some authors 
Electrically Counterpoised Dust (ECD) and has been studied with some detail, both in classical and relativistic 
theories (5lTl58j . Also, disk sources for conformastationary metrics have been considered in (25|. The method that will 
be used to construct the corresponding inverse structures is the so-called Lord Kelvin Inversion Method, applied 
to models of finite extension. The structures obtained arc extended to infinite, but one can put a cutoff due to the 
fast decay rate of the densities and their masses are finite. Moreover, the solutions obtained here agrees with all 
the energy conditions. We observe that the extremal Reissner-Nordstrom black hole can be embedded in the kind 
of spacetimes under consideration and, in particular, we note that such geometry arises naturally in the context of a 
(Kelvin) inversion as a result of a boundary term. Therefore, the superposition of a disk with the black hole turns 
out to be straightforward. This superposition can be used as a first approximation of an AGN and, due to its relative 
simplicity, may be relevant to construct more realistic models. 

The paper is organized as follows. In Sec. [ill W( 3 present the solutions of the Eisntein-Maxwell equations obtained 
in [59} . for comformastatic spacetimes. Next, in Sec. IIII1 we obtain the surface energy momentum tensor and the 
surface current density of the relativistic thin disks. In Section HVl we show the Lord Kelvin Method and we apply it 
to the solution obtained in the Sec. IIIII We verify indeed, that the solution reduces to the Minkowski one at infinity 
and we check that this is in agreement with the energy conditions. We present in Sec. IIV Al two particular families 
of solutions that represent sources of charged matter with a central hole. The distributions, which we applied this 



method to, are the Morgan and Morgan disks (MM) [7| and the flat rings introduced by [60(. Then, in Sec. IVl we 



show the way to superpose this kind of solutions with an extremal Reissncr- Nordstrom black hole and in Sec. IV Al we 
obtain the first one explicitly. Finally, in Sec IVI1 we present the conclusions of our main results. 



II. EINSTEIN-MAXWELL EQUATIONS FOR CONFORMASTATIC SPACETIMES 



For conformastatic spacetimes the line element can be written in cylindrical coordinates x M = (t, tp, r, z) as [61( 



ds 2 



e 2X dt 2 + er 2X {r 2 dtp 2 + dr z + dz 2 ), 



(1) 



where the metric function A depends on (r, z). 

The electrovacuum Einstein-Maxwell system of equations, in geometrized units such that c = G = pbQ = eo = 1, is 
given by 



T — 



F"".„ = 0, 



F F a „ F „F a P 

1 (ia J v ^Jfiu 1 up 1 - 



(2) 

(3) 
(4) 



(5) 



where is the electromagnetic four potential given by = (—0,0,0,0). Now, it is assumed that the electric 
potential <fi depends also on (r, z), in such a way that the electrovacuum Einstcin-Maxwcll system of equations 
reduces to 



V 2 A 



e- 2X V(j)- V(/», 



(0) 



2 VA • V0, 



(7) 



A,iAj = e 



-2A, 



(8) 



where i, j = 1, 2, 3, and V is the usual differential operator in cylindrical coordinates. 

A particular subclass of solutions for this system of equations was obtained in |59l | , where the authors considered a 
functional dependence between the metric function A and the electric potential 4> with respect to an auxiliary function 
U . This auxiliary function is a solution of the Laplace equation. The explicit form of A and (f> are given by 



U + k 2 ' 



(9) 



(10) 



where k\ and k 2 are integration constants, which arc chosen in such a way that the metric ((T|) reduces to Minkowski 
at infinity. As we can see when e A = 1 the electric potential is equal to zero. 



III. ENERGY-MOMENTUM TENSOR AND CURRENT DENSITY 



As we can see from ([9]) and (fT0|) . the solutions of the Einstein-Maxwell equations system ©-(O are expressed in 
terms of U, a solution of the Laplace equation. We shall deal with solutions of the Laplace equation corresponding to a 
Newtonian surface mass distribution so that, in general, U is a continuos function everywhere but its first z-dcrivativc 
must be discontinuous at z = 0. Accordingly, in order to obtain the energy-momentun tensor and the current density 
of the source, we will express the jump across the disk of the first ^-derivatives of the metric tensor as 
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and the jump across the disk of the electromagnetic field tensor as 

[F ztl ] = [A tl , z ]=2A ll , z \ z=o+ , (12) 

where the reflection symmetry of the functions with respect to z = has been used. 

Then, by using the distributional approach [62T - I641 ] . the Einstein-Maxwell equations yield the energy-momentum 
tensor 

T n» = jvvQty + T ^[i _ 0(3)] + Q^s(z), (13) 

and a current density 

J» = P l S(z), (14) 

where 9{z) and S(z) are respectively the Heaveside and Dirac distributions with support on 2 = 0. Here are the 
electromagnetic energy-momentum tensors as defined by ([3]) for the z > and z < regions respectively, whereas 

16ttQ^ = V6t - b zz 6£ + g^bl - g zz b^ + b a a {g zz 8^ - g» z 8 z v ) (15) 

gives the part of the energy-momentum tensor corresponding to the disk source, and 

47r/^ = [F MZ ] (16) 

is the contribution of the disk source to the current density. Now, the "true" surface energy-momentum tensor of the 
disk, Sf,, u , and the "true" surface current density, can be obtained through the relationships 

S,j,u = / Q^v S(z) ds n = e~ x Q^ v , (17) 



j" 



J J" 5{z) ds n = er x I^ (18) 



where ds n = ^Jg Z z dz is the proper length in the normal direction to the disk. 
Regarding the metric in Eq.([T]), the only non-zero component of Q% is 



2- 

whereas the only non-zero component of 1^ is 



p 2A A 



7° = -g. (20) 



Thus, the only non-zero component of the surface energy-momentum tensor S£ is 



S° =-^, (21) 



2tt 

and the only non-zero component of the surface current density j M is 

f = (22) 

Z7T 

where all the quantities are evaluated at z = + . 

The surface energy-momentum tensor and the surface current density of the disk, measured in rest frame of the 
dust, can be written as 

= eV»V", (23) 
f = trV», (24) 

where 

V =e- A (l, 0,0,0), (25) 
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is the velocity vector of the matter distribution. Thus, the energy density and the charge density of the distribution 
of matter are given by 

• - ^ ™ 

respectively. Now, by using Eq.([TU]), the expression ([27} can be written as 

<J = T e, (28) 

that, when static charged dust matter is included, the above conditions imply that the charge density of the disks 
must be equal to the mass density, with the possible exception of a sign. This kind of Einstein-Maxwell-extremal-dust 
solutions are generically called Majumdar-Papapetrou solutions [65], |66[ and present the interesting feature that they 
permit the construction of configurations of matter in static equilibrium. 

In particular, electrovacuum solutions with a = e = 0, in some cases, reduce to a configuration of many extremal 
Reissner-Nordstrom black holes, as was fully explored by Hartle and Hawking [6?| . Solutions with surface mass de nsity 
have been studied by (59|, and solutions with volumetric mass density have been examined in the papers of Das [5l| . 
Bonnor [52Tj55| . and others (see Ref. f56T - [58j and the references cited therein). 

Now, in the context of classical general relativity, the energy-momentum tensor is supposscd to fulfill certain 
requirements, which are embodied in the weak, strong and dominat energy conditions [f38j]. Indeed, for the case of a 
dust source, all these conditions reduce to the single condition that the energy density be greater or equal to zero, 
e > 0. On the other hand, from ^ and (f2"B")) , we have that the energy density can be written as 

where 

E=£, (30) 

is the Newtonian mass density of a disk-like source whose gravitational potential is given by U. Accordingly, if 
the Newtonian mass density S is non-negative everywhere, the corresponding rclativistic energy density e will be 
non-negative everywhere only if we take fci < 0. Therefore, in order that the energy-momentum tensor of the disks 
agrees with all the energy conditions, the Cm constants in (|33|) must be properly chosen in such a way that S > 0. 
Furthermore, if we have a Newtonian potential U that is negative everywhere, as is expected for a compact Newtonian 
source, the energy density e of the disk will be non-singular everywhere. 



IV. AXISYMMETRIC THIN DISK SOLUTIONS WITH A CENTRAL HOLE 

One of the most relevant properties of the energy density (|29[) is that it depends on r only through the potential- 
density pair corresponding to a Newtonian thin disk. So, in order to obtain solutions that correspond to rclativistic 
thin disks with a central hole we must use solutions of the Laplace equation that properly describe Newtonian sources 
with this behavior. 

A simple way to construct this kind of solutions is by using the Kelvin inversion method with finite disk-like 
systems. In the literature we can find several models representing such systems, see for example [69W751 ]. Thin disks 
of finite extension can be obtained by solving the Laplace equation in oblate spheroidal coordinates {^,rj,ip). These 
coordinates are related with the usual cylindrical coordinates by the relations 

r 2 = a 2 (l + a(l-'7 2 ), (31) 

z = <?7, (32) 
with the < £ < oo and — 1 < 7] < 1. The general solution for a disk of radius a with axial symmetry is given by [76j 

oo 

^,'?)=-E C 2nto(0^('?), (33) 
n=0 
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where P2n(v) an d 92n(£) = * 2 ™ +1 Q2n(*£) are the usual Legendre polynomials and the Legendre functions of the 
second kind respectively, and C2 n are arbitrary constants. With this general solution for the gravitational potential, 
the surface mass density takes the form 

oo 

S(r) = C ^ 2n + 1 )<hn+i(0)P2n(v)- (34) 

Now, we use the Kelvin inversion theorem that in cylindrical coordinates states that if the potential-density pair 
U (r, z), p(r, z) = Y,(r)8{z) is a solution of Poisson equation, the pair 



U(r,z) = -==U[- T — 5 ,-^— 5 ), (35) 
y/r 2 + z 2 \T + z z r l + z l ' 



p(r,z)=X(r)5(z), (36) 

where 

is also a solution of the same equation. Hence, under these transformations, the energy density given by (|29p becomes 

i(r) = k ^ = a 3 kiZ(a 2 M 

(u + h) 2 r3(fc/( ^ ;0) + fc2) 2 ' 

whereas the metric function A transforms as 

e\r, z) = J^- = y-^ , . (39) 

u + k2 7 ^u(^^)+k 2 

On the other hand, in order to have an appropriated behavior at infinity, wc must impose an additional condition 
to k\ and k 2 in such a way that we obtain an asymptotically flat spacetime. So, we will require that e A = 1 at infinity 
or 

h = fc 2 , (40) 



given that U vanishes there, as we can see in Eq. ([35]) . Moreover, under the condition in Eq. (j40|) . the electric potential 
given by (|10[) also vanishes at infinity. 

Here, if we chose solutions to the laplacc equation corresponding to finite distributions of matter without singulari- 
ties, then the Newtonian potential U is negative everywhere and the corresponding energy density will be nonsingular 
if &2 < 0. With this condition, k\ < so if the Newtonian mass distribution is always positive, it is easy to see 



that the relativistic energy density of the disk will be also positive. However, note that the factor aj \Jr 2 + z 2 
introduces singularity at the origin unless that U = at infinity. One can easily check that Eq. (|33[) implies that U 
vanishes at r/ = and £ — > oo, but in principle it is possible to add a constant to change its asymptotic behavior. 
We will assume in first place that this boundary condition is satisfied, but we shall show later that it is precisely this 
term that will give rise to an extremal Reissner-Nordstrom black hole at the center of the disk. 

Now, it is possible to show that, although the new models have infinite extension, the total masses remain finite. 
As the disks are made of dust we have that p = and the mass density is equal to the energy density. The total mass 
can be computed from the Komar integral which, in the static axisymmetric case (JT]), takes the following form 

M=2ir e(r)e- x rdr, (41) 

J a 

Now, by taking into account equations (J3SJ), f[55j) and (|40p wc find that 

(£) 3 S (V 



M = -2tt / - S ^ rdr = -2ir [ /" 2 v ' ' - rdr. (42) 

U(r) + kl J a + fci) 
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Now, wc change the variable of integration r = a / ' R, and write the mass of the disk as 

M = -2na [ = S ^ dR. (43) 

Hence, if we have a well-behaved Newtonian surface mass density and potential, e given by (|29p will be non-singular 
under the condition (|40"]) and, the total mass in Eq. (|43")) will be finite. 

In conclusion, we obtain that the general expression for the energy density, which is in complete agreement with 
all the energy conditions, can be expressed as 

*'~<ffW- (44) 

where wc have set k\ = &2 = k and with the conditions £ > and k < 0. 



A. Particular Solutions 



Now we shall restrict our general model by considering particular families of solutions by specifying the constants 
C2n in such a way that the corresponding energy densities have the characteristics we sketched previously. At first 
instance we shall deal with the MM family, which was introduced in Q (see also [zl]) and whose first member is the 
well-known Kalnajs disk [(3{j. This family are characterized by its well-behaved surface mass density, labeled with the 
positive integer m > 1, and given by 

27TCT \ a z / 

where M is the Newtonian mass of the disk and a is the radius. For each member of the family, the constants C2n 
are given by 

c = A/7r 1 / 2 (4n + l)(2m + l)! 

2n a2 2 ™+ 1 (2n+l)(m-n)!r(m + n+|) (?2n+1 (0)' v ' 



for n < m, and Cm = for n > to. So, by using this constants in the general solution (|33[) it is easy to see that the 
gravitational potential of each disk [/M will be negative everywhere, as was mentioned in the previous section. The 
closed expressions corresponding to the first four members are given by 

A(3 V 2 -1)}, (47a) 
1 OA 

-^(3r ? 2 -l) + B(357 ? 4 -30 ?? 2 + 3)], (47b) 

^(V-l) + ^W-30^+3) 

315?7 4 + 10577 2 -5)], 
20A,„ , _ 162B,_. 



u<-H>v) = 


cot x f + 

a 




cot x £ + 

a 




M r 1 , 

cot _1 ^ + 

a 




+ C(23l7? 6 - 




COt X f + 

a 








+ ^(231t? 6 




+ 69 307? 4 - 1 



(47c) 



143 



(35?y 4 - 30^ + 3) 



with 



A = 


I[(3^ + l)cot-^-3^, 






B = 


i ^[(35e 4 + 30e 2 + 3)cot- 1 C 


- 35£ 3 - 


3 ?J ' 


C = 


JL[(231£ 6 + 315£ 4 + r05£ 2 


f 5) cot" 


1 £ - 231£ 5 - 238df - 


D = 


2342912 [(6435^ + 12012^ + 


6930c: 4 - 


f 1260C 2 + 35)cot" 1 C 



35)], (47d) 

(48a) 
(48b) 

2 -f^ (48c) 



6435£ 7 - 9867c: 5 - 4213£ 3 f ] . (48d) 

35 
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It is possible to check that 



lim t/ (m) (£,0) = 0. (49) 

£—►00 



We restrict our attention to these four members but the behavior of the remaining models (to > 5) can be inferred 
from the features characterizing m = 1, 2, 3, 4. 

With the above values for the constants Cm, we can easily compute the corresponding energy density e*- 1 ™) of the 
inverted disks by using equation (|44[) . Now, in order to show graphically the behavior of each particular model, we 
first introduce dimensionlcss quantities through the relations 

t»*m = <») 

*W W . (5!) 

e (m) (r) = ae (Tn) (f), (52) 

where f = r/a, and U^ m \f) is evaluated at z — + . Accordingly, the dimensionless energy density e(r) can be written 
as 

eM(f)=- , (53) 

[£/(™)(f) + ft] 2 

with ft = aft/M. 

Then, by using the above expressions and the values of the C^n constants corresponding to the MM disks, we obtain 
the following expressions for the first four members of the family 



3ft (1 



27rf 3 



,V2 



(2r 2 - 1) - ft 



(54) 



g(2) = 



S(3) 



27T?" 3 



27rf 3 



27rf 3 



5ft (1 



.3/2 



157T 

12Nr" 



3) 



7fc(l-^ 5/2 



2 ' 



357T 



ol2 



fr (16f 6 - 24f 4 + 18f 2 - 5) - k 



9*(l-£) 



i \V2 



315ir 



32768 



(128f 8 - 256f 6 + 288f 4 - 160f 2 + 35) - k 



(55) 



(56) 



(57) 



In Fig. [TJ we plot the dimensionlcss surface energy density e^ m ^ as a function of f for the first four members of the 
family with to = 1, 2, 3 and 4. In each case, we plot e^ m ^(f) for different values of the parameter ft. The top curve of 
each plot corresponds to ft = —5 whereas the others correspond to ft = —10, —15 and —20 for the lowest curve. As we 
can see, in all these cases the energy density is everywhere positive. The curves have a central hole for f < 1, then 
it reaches a maximum and finally it decreases quickly and vanishes at infinity. Here, there are some characteristics 
of the plots that we want to point out. First, we realize that for each case, the maximum is the highest for the first 
model (m = 1) and then, it decreases as we consider the following disks (to = 2,3,...). Moreover, unlike the finite 
models studied by [59| , the behavior of the energy density with respect to the integration constant ft is the same for 
all models: e increases as |ft| decreases for all values of f, and the plots corresponding to different values of ft do not 
intersect each other. Finally, we can see that the maximum of the first model is the closest to the edge (f = 1) while 
the other curves present the maximum at a larger distance to the edge as to increases. The energy density of these 
models decays in a similar way as the Plummer-Kuzmin disks, i.e, as 1/f 3 , not too fast. The total masses can be 
computed inserting (|45jl into equation (|43]l . However, it is not possible to solve the integral analytically and hence, 
we proceeded to do it numerically for each family member. 
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FIG. 1: Dimensionless surface energy density e' m ' as a function of f for the first four members of the family with (a) m = 1; 
(b) m — 2; (c) m = 3; (d) m = 4. In each case, we plot e' m '(f) for < r < 5 with different values for the parameter k. The 
uppest curve of each plot corresponds to k = —5, and then k — — 10, —15 and —20 for the lowest curve. 




FIG. 2: In this plot, we present the mass function for the first four disk models of the family with a = 1, G = 1, k = —5 and 
M = 1. As we can see, all the masses converge to a finite value. 



In Fig. [5J we plot the total mass function for the first four models of the family. As was expected, the total masses 
are finite in spite of the fact that the models have infinite extension. 

On the other hand, we plot in Figure [3] the dimensionless surface energy density for the models m = 1, 2, 3 and 4 
keeping k at a fixed value. From this graphics it is important to notice that the models with large m decreases more 
slowly than the first models, but they have a lower maximum. Consequently we see that, in this case, each curve 
intersects the others at a certain point. 

Now let's consider disks of the same radius a and decreasing mass 

M(«) = (58) 

where S c is a constant taken equal for all disks of the MM family, n = 1,2,3,.... For this family of disks we have that 
the corresponding surface density is 



™.2 \ 2 / ^->2 



a 2 J V a 2 



n-l 



S (n) =S c 1--= . (59) 
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FIG. 3: Dimensionless surface energy density e (m) as a function of f for the first four disk models of the family with m — 1, 2, 
3 and 4. Here we plot e- m '(f) for < f < 5 setting k = — 10 for each curve. 



If we consider the superposition given by 



eW = — !!h s ( m+ i-n) ) (60) 

* — ! n (m — n) 



n=0 



we obtain the family of rings recently introduced by [6(| . We have that all these superpositions give disks of radius 
a with zero density on their centers, i.e., disks with a hole in their centers with a residual density that is smaller for 
larger m. The potentials associated to these flat rings can be found by using a superposition with the same coefficients 
as the ones used for the densities. 

For this family of models, we introduce the following dimensionless quantities 

U {m \f) = ^-ii, (62) 

£(™)(f) = (63) 

e(™)(f) = a^ m >(f), (64) 

where f = r/a, and JJ^ m ^(r) is evaluated at z — + . Accordingly, the dimensionless energy density e(f) can be written 
in the same way as equation (|53[) . but using the new definitions for U and E. 

Then, by using the above expressions for the flat rings previously introduced, we obtain the following expressions 
for the first four members of the family 
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( 8f4 + 8r 2 - 9) - fc 



(65) 



S(2) 



2 ^7 (16f 6 + 8?~ 4 + 18f 2 - 25) - fc 



(66) 



S (3) 



^-i, (640f 8 + 256f 6 + 288f 4 + 800f 2 - 1225) - fc 



(67) 



;(4) 



k\/l — Jj- 



ill 



l^rr (1792f 10 + 640f 8 + 576f 6 + 800f 4 + 2450f 2 - 3969) - fc 



(68) 





0.025 r 
0.020 - 
0.015 - 





0.5 1.0 



FIG. 4: Dimensionless surface energy density e' m ' as a function of f for the first four members of the second family with (a) 
m = 1; (b) 77i = 2; (c) m = 3; (d) m = 4. In each case, we plot e' m '(f) for < r < 3 with different values for the parameter fc. 
The uppest curve of each plot corresponds to k — —5, and then k — — 10, — 15 and —20 for the lowest curve. 



In Fig. |4l we plot the dimensionless surface energy density fS m ' as a function of f for the first four members of the 
family with m — 1, 2, 3 and 4. In each case, we plot (f) for different values of the parameter k. The top curve of 
each plot corresponds to fc = —5 whereas the others correspond to fc = —10, —15 and —20 for the lowest curve. The 
energy density is everywhere positive as was expected. The curves have a central hole for f < 1, like the first family, 
then they reach a maximum and finally they decrease very quickly and vanish at infinity. The energy density of these 
models decays as eS m ^ ~ i/f 2m+3 . As before, the total mass functions are computed numerically and in Fig. [5j we 
show the behavior of the first four models for this family. 

The behavior of these plots are very similar to those showed in Fig. [TJ They present the highest maximum for 
the model with m = 1 with the maximums decreasing as we consider the following disks. Moreover e increases as \k\ 
decreases for all values of f. Now, due to the fact that the densities decay very fast, we can put a clear cutoff and, 
thus, consider this family of models as finite flat rings. 

On the other hand, we plot in Fig. [5] the dimensionless surface energy density for the models m = 1,2,3 and 4 
keeping fc at a fixed value. In this case, we see that as m increases the maximum decreases but the variation is slower 
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Total Mass Function 




FIG. 5: In this plot, we present the mass function for the first four disk models of the family with a = 1, G = 1, k 
£ c — 1. As we can see, all the masses are finite. 



-5 and 




0.0 0.5 



FIG. 6: Dimensionless surface energy density e'"^ as a function of f for the first four disk models of the second family with 
m = 1, 2, 3 and 4. Here we plot e- m > (f) for < r < 3 setting k = —10 for each curve. 



than in the first family. Moreover, the decay is faster for larger values of m and as a result, each curve does not 
intersect the others, in contrast with the behavior of Fig. [3] 



by 



V. SOLUTIONS WITH AN EXTREMAL REISSNER-NORDSTROM BLACK HOLE 

The well-known metric of the Reissner- Nordstrom (RN) black hole spacetime in Schwarzschild coordinates is given 



{P~ P+)(P~ P-) 



dt- 



(P~P+)(P-P-) 



dp 2 + p 2 dfl 2 



(69) 



with horizons appearing at p = p± , where p± = Mhh ± \J ~ < 3 2 - ^ ne ^ me Q 2 ~ ^bh m * nc parameter space of 
electrovacuum solutions is referred to as the extremal RN spacetime. For such spacetimes the metric is reduced to 



ds 2 = _l 1 _Mu L] dt2 



1 - 



M bh V 



dp 2 + p 2 (d9 2 + sine 2 dip 2 ). 



(70) 



P J V P J 

Similar to the Schwarzschild spacetime, the RN spacetime is spherically symmetric and static. The metric (|70[) can 
be rewritten in harmonic coordinates as 



ds 2 



l + _^j dt 2 



^\ 2 [dR 2 + R 2 (d0 2 + sin 2 edip 2 )] 



(71) 



where p = R + Mbh- Some comment on the range of the radial coordinates is in order here. R > covers only 
the region outside the horizon; the region inside the horizon is obtained by continuing the solution to the range 
-M bh < R < 0. 
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Now, writing the same metric in cylindrical coordinates, we obtain 



-2 / , r \ 2 



ds 2 = -\l + // bh =) dt 2 +(l + // bh \ {dr 2 + dz 2 +r 2 d V 2 ), (72) 



\/r 2 + z 2 J V Vr 2 + 







which coincides with ([1} under the identification 



This implies that given a potential U = M^k / V r 2 + z 2 , and using Equation ([SJ, one recover the extremal RN 
spacetime. 

Now then, lets suppose that we add a constant to the potential given by ([33]) obtaining U'(r, z) = U(r, z) + C . In 
this case S'(r) = S(r) with C representing the potential at infinity. By choosing C = M^k/a, the Kelvin inversion 
theorem tell us that the inverted potential is now 

U'(r,z) = U(r,z) + -^*, (74) 

whereas the inverted density changes only by a delta factor 

S (r) = S(r) + ^ ■ (75) 
Using this new one density-potential pair we find that 

e*(r,z) = -~ (76) 

U{r, z) + k + M bh fc/Vr 2 + z 2 

and 



e-'(r) = — ^ (77) 

(C/(r,0) + fc + M bll fc/r) 2 

plus a delta term at the center. Evidently, the case U(r,z) = reduces to the metric (|72|). It is interesting to note 
that the total masses of the disks are modified by the presence of the black hole. Now M. depends on Mbh in the 
following way: 

M = -2ira [ s - - S ^ dR. (78) 

7 §U(R) + k + AI hh kaR 2 

This new solution has a singularity at the origin. From the analysis performed in section Mil it is clear that our 
solution is a non-linear superposition of an extremal RN black hole and a charged dust disk. Now, if we we look at 
the limit r — > and z — > we obtain 



or equivalcntly 



This metric has the Robinson-Bertotti form AdS2 x S 2 [78|, [79|, with anti-de Sitter (AdS) length L = Mbh, and 
coincides with the near-horizon limit of an extremal RN black hole, as expected. 
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A. An explicit solution 



In this subsection we present a fully integrated solution representing the superposition of the first member of 
inverted MM family with an extremal RN black hole. The obtention of other solutions can be done in a similar way 
and the physical properties are essentially the same. Therefore we will focus only in this model. 

Perhaps, the simplest way to obtain a closed expression for the metric coefficient is to express the potential U 
in terms of oblate spheroidal coordinates as was done in the previous sections. However, Kelvin transformation 
requires that we first express U in cylindrical coordinates which can be done analyzing the patches r < a and r > a 
separately. A more convenient expression for the potential was given recently in [75j , in which the author obtained it 
in a closed-form for all values of r in cylindrical coordinates. 

For the m = 1 disk U (r, z) is given by [7f| : 



U(r,z) 



3MG 



where 



2(2a 2 



2 , o~2\ ■ -1 / /l — h 

r + 2z ) sin 

V 2r 



+ vW/1/2 - h - *S\zWhh + h 



(81) 



/1 = yjz 2 + (r + a) 2 , 



f 2 = ^z 2 + {r-a) 2 , 

r 111 

h = a -r - z . 
On the other hand, he wave that the corresponding surface density is 



E(r) 



i^Vl-rVa 2 forr<a, 







for r > a. 



(82) 



(83) 



Therefore, the Kelvin inverted potential-density pair reduces to 



U(r,z) 



3MG 



8a 2 V~r 



where 



2 2a 2 



a\2z 2 



(r 2 + z 2 ) 2 
+ vW/1/2 - Is - 3V2 



A = 
& = 



i a 2 (z 2 + (r + 


a) 2 ) 


r 2 + z 2 


1 a 2 (z 2 + (r — 


a) 2 ) 



r 2 + z 2 
a 2 (a 2 — r 2 — z 2 ) 



{fi-&){r 2 +z 2 ) 
2a 2 r 



(84) 



(85) 



and 



S(r) 




for r < a, 
lM - y/1 - a 2 /r 2 for r > a. 



(86) 



Plugging equations and (|86p into (|76p and (|77)) . expressions for the metric coefficient and surface energy density 
of the disk can be found. 

In Figure[7J we plot the dimensionless surface energy density e as a function of f for the superposition of an extremal 
Rcissncr-Nordstrom black hole with the first member of the inverted MM family. The physical behavior of the density 
profiles are pretty similar to the plots showed in Fig. [T] In Fig. [5] we plotted the mass function for this particular 
model with different black hole masses and we found that the mass of the disk is reduced as the black hole mass is 
increased. This is precisely the expected behavior from Equation ((751) . 
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(a) (b) 




FIG. 7: Dimensionless surface energy density e as a function of r for the superposition of the first member of the inverted MM 
family with a black hole of mass (a) = M; (b) M-^h = 10M. In each case, we plot i^ m '(f) for < f < 5 with different 

values for the parameter k. The uppest curve of each plot corresponds to k = —5, and then k — —10, —15 and —20 for the 
lowest curve. 



Total Mass Function 

0.45 3- 1 1 




0.2 0.4 0.6 0.8 1 

R 



FIG. 8: In this plot, we present the mass function for the superposition of the first model with a black hole. Here we set a — 1, 
G = 1, k — —5 and M = 1 while varying Af^h- The mass of the disk is reduced as the black hole mass is increased. 



VI. CONCLUDING REMARKS 

We have presented the first superposition of a charged black hole with an annular disk made of extremal dust. To 
do this, we first obtained two infinite families of axially symmetric charged dust disks with a central hole, with well- 
behaved surface energy and charge densities. The disks were obtained by solving the electrovacuum Einstein-Maxwell 
equations system for conformastatic spacctimes, assuming a functional dependence between the metric function, the 
electric potential and an auxiliary function that was taken as a solution of the Laplace equation. Moreover, the 
solutions obtained here have a charge density that is equal, except maybe by a sign, to their mass density, in such a 
way that the electric and gravitational forces are in exact balance. 

Then, we employed the well-known Lord Kelvin Inversion Method in order to obtain annular disks, applied to the 
MM and flat rings systems, which are models of finite extension. Here is worth to mention that although this method 
was already employed by Lemos and Letelier in a previous work [3ll | , the corresponding metric functions could not be 
fully analytically computed in terms of elementary functions. The most important point to obtain the superposition 
of these disks with an extremal Rcissncr- Nordstrom black hole was the realization that such geometry arises naturally 
within the inversion by considering an extra boundary term. The solutions obtained have positive and well-behaved 
energy densities, vanishing at the edge, and as the disks are made of dust, all the models are in complete agreement 
with all the energy conditions, a fact of particular relevance in the study of relativistic thin disks models. Indeed, as 
was mentioned at the introduction, many of the relativistic thin disks models that had been studied in the literature 
do not fully agrees with these conditions. Furthermore, although the structures obtained are extended to infinity, one 
one can put a clear cutoff due to the fast decay rate of the densities, and their masses are finite. This same features 
hold for the solutions with or without the central black hole. 

We believe that the obtained models have some remarkable properties; its relative simplicity and the clear iden- 
tification of the material source associated to it not being the least of them. In particular, they may be suitable to 
study their generalizations to conformastatic models in the presence of magnetic fields as well as conformastationary 
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models with electromagnetic fields, so we are now considering some research in this direction. Also the study of the 
dynamics of test fields over these backgrounds is under consideration. 
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